matrix_product_double.cpp